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With modern scanning probe microscopes, it is possible to manipulate surface structures even at the atomic 
level. However, manipulation of nanoscale objects such as clusters is often more relevant and also more 
challenging due to the complicated interactions between the surface, cluster and apparatus. We demonstrate 
the manipulation of nanometer scale gold clusters on the NaCl(OOl) surface with a non-contact atomic force 
microscope, and show that the movement of clusters is in certain cases constrained to specific 
crystallographic directions. First principles calculations explain this kinetic anisotropy as the result of the 
cluster attaching to surface defects: cation vacancies allow the clusters to bond in such a way that they only 
move in one direction. Constraining the movement of clusters could be exploited in the construction of 
nanostructures or nanomechanical devices, and the manipulation signatures may also be used for 
identifying cluster-defect complexes. 
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Controlling structures at the atomic level is the ultimate challenge in nanofabrication - absolute control over 
the positions of individual atoms, molecules, or clusters on substrates allows for nearly complete freedom 
in material composition or in the fabrication of nanodevices. Although arranging structures atom-by- 
atom is attractive from the perspective of fundamental interest in surface physics and chemistry, for many 
applications manipulation of larger units is more relevant. With respect to metallic clusters, interest stems from 
catalysis where clusters promote many catalytic reactions 1 " 3 , but also from nanotechnology where clusters can be 
used as building blocks for nanoelectronic devices 4 . In tribology, clusters can even serve as units for probing the 
friction at the nanoscale 5 ' 6 . In all these domains, manipulating clusters as single objects enables the systematic 
study of the impact of the specific and voluntarily chosen adsorption site on the physical, mechanical, electronic 
and chemical properties of a single cluster 7 . While manipulation of atoms was achieved already in 1990 with 
scanning tunneling microscope (STM) 8 , atomic manipulation with non-contact atomic force microscopy (nc- 
AFM) 9 has been a much more recent development 10 . In particular, manipulation on insulating surfaces remains 
particularly challenging 1112 , yet insulators are a crucial class of materials in nanoelectronics and nanocatalysis. 

Similar to STM 13 , two methods are used to manipulate clusters by nc-AFM: the first method consists of moving 
the cluster by pushing during a spectroscopy step (Fig. 1 (a)). Here, the AFM tip is brought down vertically near 
the cluster until at some point the repulsive interaction between the tip and the cluster is strong enough to make 
the cluster move. In the second method the tip is scanned above the cluster (Fig. 1 (b)), preferentially in the 
constant height AFM mode 14 " 16 . At a specific distance, as the tip approaches the cluster laterally, a force between 
both moves the cluster. For simplicity, we will call the former method kicking and the latter sliding. 

This work investigates methods for directing the manipulation of a nanocluster beyond the constraint of the 
tip's scanning direction. In particular we demonstrate the possibility of using intrinsic defects on the surface to 
mediate the diffusive behaviour of the clusters. We have chosen the NaCl(OOl) surface as a model insulator 
substrate firstly because it can be easily prepared and examined by nc-AFM even at the atomic scale 17 , but also 
because the characteristic defects at it's surface are well-established 18 . Furthermore, NaCl and other alkali halide 
surfaces have been key in breakthrough nc-AFM studies in atomic 19 " 21 and charge identification 18 ' 22 , as well as 
providing prototype substrates for the study of adsorbed molecules 23 " 25 , self-assembly 26 , friction 27 ' 28 , and local 
electronic structure 29 ' 30 . 

Results 

Manipulation on terraces. When single neutral gold atoms are deposited on perfect NaCl(OOl) surfaces, the gold 
atoms diffuse on the terraces due to their high mobility at room temperature 31 " 33 until they reach a step-edge, 
where they aggregate into clusters 1 ' 34 . At room temperature and low gold coverage, the steps are decorated by gold 
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Figure 1 | Principles used in the manipulation of clusters, (a) the 

cluster is moved by kicking during a spectroscopy step; (b) moving the 
cluster by sliding during scanning. 



clusters whereas on the flat terraces very few clusters can be found 
because of the lack of surface defects that would act as nucleation sites 
(see further discussion in the Supplementary information). The 
image in Fig. 2 (a) exemplifies a typical decoration of steps (along 
(100) equivalent directions) by gold clusters. The clusters have a 
mean diameter of 5 nm (—2500 atoms), which was estimated 
from the flux of incoming gold during the deposition. 

In this first series of manipulation experiments, we focus on one of 
the rare clusters on the terrace (highlighted in green in Fig. 2) and 
consider its manipulation characteristics. The cluster was first 
imaged with a high magnification in the constant height mode 
(Fig. 2 (b)). In the figures, the slow scanning direction of the AFM 
tip is always from the bottom to the top. The same cluster was imaged 
several times while gradually bringing the tip closer to the surface. At 
a specific tip-surface distance the cluster suddenly disappeared dur- 
ing imaging (arrow in Fig. 2 (e)) by escaping behind the tip to the 
region that had just been imaged to the position shown in Fig. 2 (c). 
This kind of backheel led to a 16 nm long movement of the cluster 
along the [100] direction. In a second manipulation step the same 
cluster was imaged again and the tip -surface distance was further 
reduced. Now, the cluster could be moved by sliding in front of the tip 
along the [110] direction (Fig. 2 (f)). Because the cluster is moved in 
front of the tip from line to line, it appears with an elongated contrast, 
reflecting its movement. The new position of the cluster, which is 
11 nm further away from the old one, can be seen in Fig. 2 (d). 

Manipulation from step-edges. In Fig. 3, similar clusters of about 
2000 atoms decorate intersecting (110) and (100) steps on the 
NaCl(OOl) surface. One of the clusters is highlighted (I in Fig. 3 
(a), green color) at the (110) step-edge. After recording a series of 
high resolution images at this site, the last image of this series shows 
that the cluster was moved 12 ± 1 nm from the step (to position II in 
Fig. 3 (b)). We continued to move the same cluster, now on the 
terrace, by sliding about 15 ± 1 nm along the (110) direction 
(Fig. 3 (c)), which has been verified in a following image with a 
larger scale (Fig. 3 (d)). We could move the cluster also by kicking 
with a vertically approaching tip in a spectroscopy step. Figure 3 (e) 



shows that in this case the cluster could also be moved by 16 ± 1 nm 
from its initial position. Altogether, the cluster was moved by the tip 
dozens of times during the whole manipulation process. Regardless 
of the method of manipulation, the events always resulted in cluster 
movement along equivalent (110) directions on the surface (see path 
of cluster in Fig. 3 (f)). This extremely anisotropic movement is in 
clear contrast to the cluster found on the terrace, which could easily 
be moved in different directions. 

Calculated manipulation kinetics. In order to understand the 
movement of clusters during manipulation and pinpoint the origin 
of the experimentally observed differences, we sample the possible 
pathways for moving 20 and 100 atom gold clusters (Fig. 4 (a)) using 
density functional calculations. With the smaller Au 20 cluster we 
perform extensive analysis on the adsorption configurations and 
diffusion pathways while the larger Au 100 cluster is used in selected 
cases to assess scaling with respect to cluster size. 

The Au 2 o cluster is found to adsorb on the perfect surface (adsorp- 
tion energy 0.5 eV) and the energy barriers for moving the cluster in 
any given direction are very low, about 0.2 eV (Fig. 4 (b)). At room 
temperature, such a diffusion barrier corresponds to a Boltzmann 
factor of 10" 4 , and thus such a small cluster should move sponta- 
neously due to thermal energy. Furthermore, the barriers for moving 
the cluster in either (100) or (1 10) directions are practically the same 
and thus no direction should be preferred. 

The calculated high mobility on the ideal terrace strongly implies 
the presence of a defect under the cluster in Fig. 2 — else the cluster 
would have diffused to a step during growth. Hence, we also consider 
the role of defects in cluster adsorption and diffusion. The calcula- 
tions show that in the presence of surface vacancies, or F-centres, the 
cluster adsorbs even more strongly than on the perfect surface, as the 
adsorption energies are 2.1, 2.5, and 1.1 eV for adsorption at neutral 
CI, Na and double vacancies, respectively. (Defect charge state influ- 
ence is discussed in the Supplementary information.) Also the dif- 
fusion energy barriers are much higher when defects are present, as 
discussed below, meaning it is much more difficult to move the 
cluster, and this verifies that defects can act as anchors for the clusters. 

In the presence of a CI vacancy, the barriers for moving the cluster 
are 0.8 - 1.1 eV and the anisotropy is relatively weak (Fig. 4 (c)). A 
cluster bonded to a Na vacancy, on the other hand, prefers a [110] 
direction where the barrier is calculated to be 0.7 eV. In comparison, 
the (100) barriers are 0.9 eV, and in the perpendicular [1 10] dir- 
ection the barrier is 1.2 eV (Fig. 4 (d)). Similar calculations with a 
100 atom gold cluster show that this is still the case even when the 
clusters grow. Figure 4 (e) shows the calculated barriers for the Au 10 o 
+ Na vacancy complex. When the vacancy is located near the cluster 
edge, in the energetically favored position, the anisotropy is even 
stronger than for the small cluster (barriers 0.5 - 1.1 eV). 



step »» % ^ 
un V^tllO] 

Figure 2 | Manipulation of a Au cluster (green) on NaCl(OOl) along the [110] and [100] directions (nominal Au thickness: 0.3 mono-layers, 

flux: 4.1 X 10 12 atoms/cm 2 /s). (a) The surface region before the manipulation (density: 0.9 X 10 11 clusters/cm 2 , cluster- cluster mean distance at steps: 

~ 20 nm). (b) - (d) The cluster before (b), after the first (c) and after the second (d) manipulation experiment, (e) During the first manipulation, a part of 

the cluster could be still imaged until the position labeled by the arrow where the cluster suddenly started to escape behind the tip (backheel). (f) During 

the second manipulation step, the cluster was moved by sliding along the [110] direction, (g) Schematic of the two manipulation steps. 

(a) 250 X 250 nm 2 ,A/= -22 Hz, (b,c,d) 65 X 136 nm 2 , A/« -22; -18; -20 Hz, (e) 40 X 40 nm 2 , Af= -48 Hz, (f) 30 X 26 nm 2 , Af= -69 Hz, All 

images: f 0 = 275 kHz, v Scan = 10 Hz (e: 24 Hz, f: 33 Hz), k = 28 N/m, A pp = 15 nm. 
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Figure 3 | Experimental constant height images presenting several manipulation steps of Au clusters on NaCl(OOl) (nominal thickness: 0.64 mono- 
layers, flux: 1.0 X 10 13 atoms/cm 2 /s). (a) Decoration of NaCl steps by Au clusters (cluster-cluster distance: 8.4 nm, cluster density: 4.6 X 10 11 
clusters/nm 2 ). (b) A single cluster (green) is separated from the step during a scan with the AFM. The dashed circle represents the position of the 
manipulated cluster in the previous image, (c) Controlled manipulation of the cluster by sliding, (d) The result after the movement, (e) Controlled 
manipulation of the cluster by kicking in a spectroscopy step, (f) Schematic of the manipulation steps, (a) 88 X 88 nm 2 , Af= — 52 Hz, 
(b, d, e) 54 X 51 nm 2 , A/« -54 Hz, (c) 20 X 20 nm 2 , Af= - 120 Hz, All images: f 0 = 247 kHz, k = 32 N/m, v Scan « 15 Hz, A pp = 14 nm. 



Discussion 

The charged Na + and CI" vacancies (latter are F + -centres) are 
expected to be mobile on the NaCl(OOl) surface and diffuse to steps 
during annealing 35 ' 36 . Recent nc-AFM and Kelvin Probe Force 
Microscopy (KPFM) studies have further shown that Na + vacancies 
are the dominant defect at step -edges and kinks 18,37 , which is due to 
the Debye-Frenkel layer formation on alkali halide (001) surfaces 38,39 . 
This suggests that a cluster taken from a step, as shown in Fig. 3, is 
with a high probability anchored by a Na + vacancy. The experiment- 
ally observed kinetic anisotropy of these clusters can then be 
explained as the characteristic behaviour of this particular cluster- 
defect complex, as revealed by the simulations. The third type of 
possible defect on NaCl(OOl) is the F°-center (CI" vacancy + elec- 
tron, i.e., a Cl° vacancy), and these defects are expected to be stable on 
the terrace (calculated diffusion barrier for the defect is 0.9 eV). 
Therefore, they are also the most probable defect on which clusters 
found on terraces are adsorbed (as in Fig. 2). From this consideration, 
the theoretically predicted small diffusion anisotropy for a cluster 
adsorbed on a CI vacancy agrees with the experimentally found 
isotropic movement of the cluster on the terrace. 

The reason why the defects influence the movement of clusters 
differently can be understood by examining the calculated structures 
of the bonded cluster- vacancy pairs. In Fig. 5 (a) and (b), the differ- 
ential charge densities are shown for clusters in the presence of CI 
and Na vacancies, respectively. This density represents the change in 
electronic charge upon bringing the cluster and the surface 
together 40 . In the presence of a CI vacancy, the Au cluster bonds 
almost exclusively at one of its corners, attaching to a CI ion next 
to the vacancy. In fact, the cluster tilts 10° as the bonded corner is 
attracted to the CI ion. Because of this, the cluster prefers to move by 
pivoting around the bonded corner, as schematically shown in Fig. 5 
(c) (although the preference is not strong, as shown by the almost 
isotropic barriers in Fig. 4 (c)). Simultaneously, the CI vacancy jumps 
to a neighboring site, allowing the cluster to bond to another corner. 
Combining these rotations leads to movement in the (100) direction, 
but since the cluster is constantly reorienting itself, on average none 
of the (100) surface directions is preferred and the cluster may in fact 
move in any direction by taking a zig-zag path. 

On the other hand, if the cluster attaches to a Na vacancy, it does 
not bond at a corner, but along an edge which orients itself with an 
entire row of CI ions on the surface. This is seen in Fig. 5 (b). Again, 
the cluster is seen to slightly tilt (5.1°) towards the active edge to 
facilitate the bonding. For such a configuration, the energetically 
favored mode of movement of the cluster is sliding along the row 
of CI ions with which it has bonded, so that the active edge always 



stays on the CI row. This is represented in Fig. 5 (d). Since this 
movement follows a single row of CI ions, the cluster moves in a 
[110] direction. Notably, since the cluster cannot move perpendicu- 
larly to the CI row without detaching or rotating first, a much more 
difficult task than sliding in the preferred direction, this configura- 
tion results in truly anisotropic movement. Exactly the same kind of 
attaching to a CI row occurs with a larger cluster as shown in Fig. 5 (f), 
demonstrating that this behaviour does not depend on the cluster 
size. 

There is a discrepancy between the cluster sizes in the presented 
theory and experiments (100 and 1000 atom clusters, respectively). 
Although such a difference in the number of atoms only means a 
linear size difference by a factor of two, difference in size at this scale 
may lead to changes in the cluster morphology. On ionic surfaces like 
NaCl(OOl) or MgO(OOl), gold clusters expose either their (001) or 
(111) facets to form rectangular or triangular shapes, respectively, 
and the shape may depend on the size and growth conditions 44,45 . 
However, as our simulations predict the anisotropic movement for 
the cluster- vacancy complex to be due to bonding of low-coordi- 
nated gold atoms at the cluster edge with CI ions on the substrate near 
a Na vacancy, the anisotropy should not depend on the size and exact 
shape of the cluster. Instead, the behaviour is a characteristic feature 
of this kind of cluster- defect pair. We therefore expect the same 
characteristics also for the larger clusters, as seen in the experiments. 

To summarize, we experimentally demonstrate nanomanipula- 
tion of gold clusters roughly of the size of 2000 atoms (5 nm dia- 
meter) as single units on the terraces of the insulating NaCl(OOl) 
surface using nc-AFM. From first principles calculations we show 
that stable clusters on the terrace are expected to be adsorbed at 
defects, which always accompany a cluster during its movement on 
the surface. Since the NaCl(OOl) surface only has vacancies as poten- 
tial defect candidates at the surface, the manipulation properties are 
dominated by the type of vacancy: clusters on cation vacancies 
exhibit an anisotropic movement along the (110) surface directions 
whereas clusters on anion vacancies do not. Especially with respect to 
the cation vacancies, which are known to be located at the step edges, 
theory perfectly agrees with experiments since clusters that are 
detached from steps by the tip move indeed along (110) surface 
directions. 

In general, this work demonstrates both experimentally and the- 
oretically that even on smooth surfaces cluster- defect pairs can have 
extremely anisotropic movement patterns, and that this anisotropy 
can be strong enough to be the dominant factor controlling their 
manipulation. Such constrained movement can be exploited when 
building structures through nanomanipulation or even in the design 
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Figure 4 | Calculated diffusion barriers for various cluster- vacancy 
structures in eV. (a) The Au 2 o and Au 10 o clusters used in calculations, 
(b) Au 2 o without defect (c) Au 2 o on a CI vacancy (d) Au 2 o on a Na vacancy 

(e) Au 10 o on a Na vacancy (vacancy at the energetically favoured edge site). 

(f) The (110) directions are set to run horizontally and vertically in the 
figures. Positions of Na and CI vacancies ( V Na and V cl ) are denoted by dark 
and light squares, respectively. 

of nanomechanical machines. Also conversely, observing specific 
anisotropy patterns could be used for distinguishing the types of 
defects on which clusters are adsorbed. Of course, this is not neces- 
sarily a one-to-one mapping: in a system where different kinds of 
defects exist, it is possible that different cluster- defect complexes have 
very similar manipulation patterns. Still, examining the cluster move- 
ment may allow one to distinguish the most likely defect types, espe- 
cially if the intrinsic occurrence probabilities of defects are known. 

Aside from providing an important characterization of defects, the 
ability to control and create defect- cluster complexes offers great 
potential for future studies in catalysis where, for instance, the cluster 
reactivity must be studied as a function of adsorption site. Such 
manipulation experiments can be combined with a control over 
the defects, which can be created by e.g. the AFM tip followed by 
imaging and spectroscopy identification, or by doping or irradiation 
of the substrate material such that known defects appear on the 
surface. When manipulation can be used for an overall characteriza- 
tion of the defect- cluster system, it will be a powerful tool for linking 
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Figure 5 | Differential charge densities /^u+Naci — Au — /faaci plotted at 
isosurfaces —0.01 eA~ 3 (grey) and 0.01 eA~ 3 (orange) for Au 20 on a CI 
vacancy (a) and on a Na vacancy (b). The [110] direction runs into the 
image. Schematic representation of the preferred diffusion modes for the 
Au 20 on a CI vacancy (c) and a Na vacancy (d), with the [110] direction of 
the NaCl running vertically as shown in (e). Differential charge density for 
Au 10 o on a Na vacancy (f) at the same isosurfaces as in (a) and (b). 
Rendered using VMD 57 . 



observed properties with the atomic and electronic structure of the 
defect- cluster complex. 

Methods 

Sample preparation. The undoped NaCl crystals are grown in an argon atmosphere 
by the Czochralski growth method in our laboratory at the CINaM 41 . Before surface 
preparation, crystals and sample holders are cleaned by annealing (~250°C) in an 
oven located inside the main UHV chamber 42 . After cooling, the (001) surfaces are 
prepared by cleavage of the crystals along the (001) plane at room temperature in the 
UHV chamber 42 . The crystals are then annealed a second time in the same oven at 
~ 150°C for a few hours to remove surface charges 43 . The preparation of a clean 
NaCl(OOl) surface by cleavage and following annealing, the deposition of gold onto 
such a NaCl(OOl) surface, and following AFM experiments are done in one UHV 
chamber, which is equipped with a sample transfer system. 
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Gold coverages were estimated by measuring the flux of a beam of neutral gold 
atoms from a Knudsen cell using a quartz crystal microbalance. The pressure during 
the deposition is in the lower 10~ 9 mbar range. The nc-AFM images were not used for 
coverage estimation due to the strong tip-cluster convolution, which is of importance 
when the clusters have a size of only a few nanometers 1516 . Especially for the small 
clusters studied in this work, their size (diameter ~ 5 nm) is smaller than the tip apex, 
the diameter of which is typically < 20 nm for new unused tips 46 . In this case, the 
shape of the tip apex is basically imaged at each cluster site, which becomes imme- 
diately evident when the apex changes upon tip changes during imaging. The con- 
volution therefore overestimates the cluster size and induces large errors in any kind 
of estimates where the cluster size is a parameter. 

AFM and manipulation protocols. Frequency- modulated noncontact AFM (nc- 
AFM) and Kelvin probe force microscopy (KPFM) experiments are performed in the 
low 10" 10 mbar pressure range and at room temperature with an Omicron RT-AFM/ 
STM. Conducting silicon cantilevers are used (Nanosensors, p-Si, 0.015 Q.cm, 247 
and 275 kHz resonance frequency, 32 (247 kHz tip) and 28 N/m (275 kHz tip) 
spring constant). Noncontact AFM primarily measures the change of frequency 
(detuning Af) of the cantilever oscillation that is due to the interaction of the attached 
sensing tip with the surface 17 ' 47 . In standard constant Af operation mode, a high 
feedback gain of the distance regulation loop is used such that variations in Af from a 
pre-set value ( A/ pre . set ) are nullified by regulating the tip-surface distance (topography 
images). For imaging the surface in a quasi-constant height mode, which was the 
scanning mode used for this work, the gain is drastically lowered such that the tip 
follows only the mean inclination of the surface. In this mode, the tip-surface 
interaction is stored in the Af image. When scanning above clusters the tip is closest 
above the top facet of the clusters so that the tip -surface interaction is strongest (dark 
contrast = stronger interaction = more negative detuning Af). More details about 
scanning modes with respect to clusters can be found in Refs. 15, 16. In order to 
measure the detuning, a demodulator from NanoSurf is used (EasyPLL). All images 
presented in this work were conducted in the constant height mode of the AFM. 
Images were acquired with the Omicron SCALA system and analyzed with the 
Gwyddion software 48 . 

Two methods are used to manipulate the clusters by nc-AFM. The first method 
consists of moving the cluster by pushing during a spectroscopy step (Fig. 1 (a)). The 
scanning is switched off, the tip is positioned in a specific spot located on or very close 
to the cluster and the tip-sample distance is finally decreased monotonously. At a 
specific tip-cluster distance, at which the tip apex is very close to the cluster, the cluster 
moves on the surface due to a sufficiently large lateral force between the foremost 
atoms of the tip and the cluster. In the second method the tip is scanned preferentially 
in the constant height mode 14 " 16 above the cluster (Fig. 1 (b)). From image to image, 
the tip-surface distance is reduced by choosing more negative detuning values Af. At a 
specific distance, the tip comes from the side and is sufficiently close to the cluster 
such that again a lateral force between both moves the cluster. This produces an 
elongated contrast of the cluster, as can be seen in Fig. 2(f) and Fig. 3(c). From this 
specific contrast it can be concluded that the cluster slides by some small distance in 
front of the tip - the tip images almost always the same part of the cluster such that the 
same contrast is reproduced in each scanning line. 

Calculations. Calculations were performed using spin polarized density functional 
theory (DFT), as implemented in the periodic plane wave code vasp 49 ' 50 . Generalized 
gradient functionals (PBE) 51 and the projector augmented wave method (PAW) 49 ' 52 
were applied. An energy cutoff of 300 eV and a 2 X 2 X 1 /c-point mesh were found to 
converge all key parameters of the system. The NaCl (001) surface was modeled using 
slabs three ionic layers thick with the bottom layer frozen, separated by 2.5 nm of 
vacuum. Energy barriers were determined using the climbing image nudged elastic 
band method 53 , where typically three or five chain images were sufficient. For 
simplicity, we ignore the effect of the AFM tip and concentrate on examining the 
energy landscape as felt by the clusters on the NaCl (001) surface. 

Dispersion forces were not explicitly included in the calculations, and although 
they may influence adsorption energies, their effect on the diffusion barriers is 
expected to be small 54 ' 55 . In order to verify that this is the case, we calculated as a 
benchmark a series of diffusion pathways using both PBE and the vdW-DF func- 
tionals 56 , the latter implementing van der Waals forces. We especially used the vdW- 
DF functional to re-examine the diffusion pathways predicted by PBE to have the 
lowest and highest barrier energies for the Au 20 cluster and the different defects. In all 
cases, the functional including dispersion forces yielded energy barriers close to that 
given by PBE, 0.0 - 0.2 eV lower, and the hierarchy of energies was the same for both 
functionals. 

Lateral forces needed to move the clusters were also estimated from the calculated 
trajectories. In general, in the preferred pathways where the energy barriers were 
found to be the lowest, also the lateral force on the cluster was low, in the range of 
0.1 - 0.2 eV/A. On the high barrier trajectories also the lateral forces were higher, 
about 0.2 - 0.3 eV/A. That is, a direct correlation was seen between the energy 
barriers and lateral forces. We focus on the energy barriers in the discussion, since 
these are more accurately determined in the calculations. 

The calculated results were obtained for a charge neutral system unless stated 
otherwise. This means that to create a Na vacancy, we remove a Na atom from the 
lattice leaving behind effectively a Na + vacancy and a missing electron. Similarly, a CI 
vacancy is effectively a Cl~ vacancy with an additional electron. In the case of the 
anion vacancy, an electron remains localized in the vacancy forming an F°-centre. 
This is the expected form of defects at the surface, particularly with the abundance of 



electron sources (i.e. metal) in the system. However, it is also possible to change the 
number of electrons in the calculations to simulate Na + and CI" vacancies. These 
systems are not charge neutral, but the excess is compensated by a background charge. 
Extensive tests were performed to study the influence of the charge state of the defects 
on the results. These calculations are discussed in the Supplementary information. 

20 and 100 atom gold clusters were used in the simulations, constructed in the 
shape of a tetrahedral pyramid or a truncated pyramid, respectively, as shown in Fig. 4 
(a). Much larger systems are not feasible at this level of accuracy. Gold clusters prefer 
to grow in (100) epitaxy but also in (1 1 1) epitaxy with a triangular or hexagonal 
shape 1 ' 44 ' 45 . For consistency, we focus on small Au 20 and Auioo clusters in (111) 
epitaxy in the calculations. For Au 100 > our calculations show the (111) faceted trun- 
cated tetrahedron structure to be energetically preferred to the best (100) epitaxial 
geometry we could find by 2.5 eV. 
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